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CN ; ABSTRACT 

In the present work, by the help of the newly released Union2 compilation which consists of 557 
Type la supernovae (SNIa), we calibrate 109 long Gamma-Ray Bursts (GRBs) with the well-known 
' Amati relation, using the cosmology-independent calibration method proposed by Liang et al.. We 

have obtained 59 calibrated high-redshift GRBs which can be used to constrain cosmological models 
without the circularity problem (we call them "Hymnium" GRBs sample for convenience). Then, we 
consider the joint constraints on 7 cosmological models from the latest observational data, namely, 
the combination of 557 Union2 SNIa dataset, 59 calibrated Hymnium GRBs dataset (obtained in 
this work), the shift parameter R from the WMAP 7- year data, and the distance parameter A of the 
measurement of the baryon acoustic oscillation (BAO) peak in the distribution of SDSS luminous 
red galaxies. We also briefly consider the comparison of these 7 cosmological models. 
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I. INTRODUCTION 

Since the discovery of current accelerated expansion of our universe Type la supernovae (SNIa) 
have been considered to be a powerful probe to study this mysterious phenomenon. However, SNIa are 
plagued with extinction from the interstellar medium, and hence the current maximum redshift of SNIa is 
only about z ~ 1.755. On the other hand, the redshift of the last scattering surface of cosmic microwave 
background (CMB) is about z ~ 1090. There is a wide "desert" between the redshifts of SNIa and CMB. 
So, the observations at intermediate redshift are important to distinguish cosmological models. 

Recently, Gamma-Ray Bursts (GRBs) have been proposed to be a complementary probe to SNIa (see 
e.g. Q and references therein). So far, GRBs are the most intense explosions observed in our universe. 
Their high energy photons in the gamma-ray band are almost immune to dust extinction, in contrast to 
supernovae. Up to now, there are many GRBs observed at 0.1 < z < 8.1, whereas the maximum redshift 
of GRBs is expected to be 10 or even larger [3j. Therefore, GRBs are considered to be a hopeful probe to 
fill the "desert" between the redshifts of SNIa and CMB. We refer to e.g. 0,0^1 for some comprehensive 
reviews on the so-called GRB cosmology. As is well known, there is a circularity problem in the direct 
use of GRBs 0, mainly due to the lack of a set of low-redshift GRBs at z < 0.1 which are cosmology- 
independent. To alleviate the circularity problem, some statistical methods have beenproposed, such as 
the scatter method @ , the luminosity distance method @ , and the Bayesian method Q . Other methods 
trying to avoid the circularity problem have been proposed in e.g. [id lllj . Recently, a new idea of the 
distance ladder to calibrate GRBs in a completely cosmology-independent manner has been proposed 
in [TH, HH independently. Similar to the case of calibrating SNIa as secondary standard candles by using 
Cepheid variables which are primary standard candles, we can also calibrate GRBs as standard candles 
with a large amount of SNIa. And then, the calibrated GRBs can be used to constrain cosmological 
models without the circularity problem. We refer to e.g. [IH, [lj, [l5[ for some relevant works. 

It is worth noting that in the literature most relevant works mainly used the 69 GRBs compiled in Q 
or the 70 GRBs compiled in [l6j]. As of the end of 2009, there are 109 long GRBs with measured redshift 
and spectral peak energy compiled in [l7j. The number of usable GRBs is significantly increased (60j . 
Among these 109 GRBs, the data of 70 GRBs are taken from [iff; the data of 25 GRBs are taken 
from 0; the data of remaining 14 GRBs (090516, 090618, 090715B, 090812, 090902B, 090926, 090926B, 
091003, 091018, 091020, 091024, 091029, 091127, 091208B) are provided by L. Amati in private com- 
munication On the other hand, the cosmology-independent method to calibrate GRBs proposed 
in [l2j needs a large amount of SNIa at z < 1.4. Thus, the SNIa dataset plays an important role in the 
calibration of GRBs. Previously, the relevant works (e.g. [H, [HI, [l5| ) used the 192 Davis SNIa dataset or 
the 307 Union SNIa dataset to calibrate GRBs. Very recently, the Supernova Cosmology Project (SCP) 
collaboration released their Union2 compilation which consists of 557 SNIa [20| . The Union2 compilation 
is the largest published and spectroscopically confirmed SNIa sample to date. The number of usable 
SNIa is also significantly increased. Obviously, a larger SNIa dataset could bring a better interpolation 
for the cosmology- independent method to calibrate GRBs proposed in [l2j|. Motivated by both the sig- 
nificant improvements in GRBs and SNIa, in the present work, we update the calibration of GRBs with 
the well-known Amati relation, using the cosmology-independent method proposed in [l2l | , and obtain 59 
calibrated high-redshift GRBs which can be used to constrain cosmological models without the circularity 
problem (we call them "Hymnium" GRBs sample for convenience). 

On the other hand, very recently the WMAP Collaboration also released their 7-year CMB data 



(WMAP7) in [21|. Therefore, it is natural to consider the joint constraints on cosmological models from 
these latest observational data, namely, the combination of 557 Union2 SNIa dataset [20|, 59 calibrated 
Hymnium GRBs dataset (obtained in this work), the shift parameter R from the WMAP7 data [2l|, and 
the distance parameter A of the measurement of the BAO peak in the distribution of SDSS luminous red 
galaxies [22, 23j. In this work, we consider 7 cosmological models and obtain the observational constraints 
on them. Notice that we assume the universe to be spatially flat throughout this work. 

This paper is organized as followings. In Sec. [Hi by the help of the newly released Union2 compila- 
tion which consists of 557 SNIa, we calibrate the 109 GRBs with Amati relation, using the cosmology- 
independent method proposed in [l2| . In Sec. IIII1 we briefly introduce the methodology to constrain 
cosmological models with the combined latest observational data. Then, we obtain the corresponding 
constraints on 7 cosmological models in Sec. IIVI Further, we consider the comparison of these 7 models 



in Sec. [V] Finally, we give the brief conclusions in Sec. [VI 
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II. UPDATING THE CALIBRATION OF GRBS WITH UNION2 SNIA DATA 

In this section, as mentioned above, we are going to calibrate the 109 GRBs [13 (and 0, El, El) 
by using the 557 Union2 SNIa data poj . In fact, we are closely following the cosmology-independent 
calibration method used in [l2|, Ell- Notice that in the present work, we calibrate GRBs only with 
the Amati relation, so that we can use a larger GRBs dataset for single GRB luminosity relation. As 
in [IHEll, we choose z = 1.4 to be the divide line to separate GRBs into low- and high-redshift groups. 
In the 109 GRBs compiled by Amati [13 (and [MEiGjl), there are 50 GRBs at z < 1.4 and 59 GRBs 
at z > 1.4. The maximum redshift of these 109 GRBs is z = 8.1 for GRB 090423. 

Several years ago, Amati et al. found the well-known E p ^ — Ei SO correlation in GRBs as E Pj i = K X E™ Q 
by using 12 GRBs with known redshifts 24] (note that this correlation has been firstly discovered by 
Lloyd et al. [62j independently), where E Pt \ — E Pi0 b s x (1 + z) is the cosmological rest-frame spectral 
peak energy; the isotropic-equivalent radiated energy is given by 

E iso = 4nd 2 L S holo (l + z)- 1 , (1) 

in which Sboio is the bolometric fluence of gamma rays in the GRB at redshift z, and is the luminosity 
distance of the GRB. Later, Amati et al. have updated it in [T614T8I . |25|. Up to now, some theoretical 
interpretations have been proposed for the Amati relation It might be geometrical effects due to 
the jet viewing angle with respect to a ring-shaped emission region J26[, or with respect to a multiple 
sub-jet model structure [13 • An alternative explanation of the Amati relation is related to the dissipative 
mechanism responsible for the prompt emission [281 ] . However, we alert the readers to the fact that there 
is a considerable debate as to whether the Amati relation is an intrinsic effect or the result of detection 
biases (or a combination of these two) (63[. We refer to e.g. [6l|, [64j for some notable works which discuss 
the Amati relation in particular. They noted that many outliers are not included in the studies of Amati 
and Schaefer, and argued that this relation might simply be wrong, or its inferred scatter might be too 
low or different from the observed scatter [_63j - Anyway, let us go ahead. For convenience, similar to Q, 
we can rewrite the Amati relation as 

log^ = A + 61og-^i-, (2) 
erg 300 kev 

where "log" indicates the logarithm to base 10, whereas A and b are constants to be determined. In 
the literature, the Amati relation was calibrated with Ei so computed by assuming a ACDM cosmology 
with specified model parameters. As mentioned above, this is cosmology-dependent and the circularity 
problem follows. Here, we instead use the method proposed in 1121 to calibrate the Amati relation in a 
cosmology-independent manner. We refer to the original paper 12] for more details. 

Firstly, we derive the distance moduli for the 50 low-redshift (z < 1.4) GRBs of [13 (and [11 [II 01]) 
by using cubic interpolation from the 557 Union2 SNIa compiled in [20j. We present the interpolated 
distance moduli /i of these 50 GRBs in the left panel of Fig. [T] The corresponding error bars are also 
plotted. In TablelH we give the numerical data of these 50 GRBs. As in [l2j], when the cubic interpolation 
is used, the error of the distance modulus /i for the GRB at redshift z can be calculated by 

where A { = Y[ {zj _ z ) / H ( z j ~ z i) > 

in which j runs from 1 to 4 but j ^ i; on the other hand, e^j are the errors of the nearby SNIa whose 
redshifts are Zi. Then, using the well-known relation 

, = 51og^ + 25, (3) 

one can convert distance modulus n into luminosity distance (in units of Mpc). From Eq. ([T]) with the 
corresponding Sboio given in [H, H; EH j we can derive E[ so for these 50 GRBs at z < 1.4. We present 
them in the right panel of Fig. [I] whereas E Pt i for these 50 GRBs at z < 1.4 are taken from [l6|, El, EH- 
Also, we present the errors for these 50 GRBs at z < 1.4, by simply using the error propagation. From 
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Fig.[T one can clearly see that the intrinsic scatter is dominating over the measurement errors. Therefore, 
as in 1, the bisector of the two ordinary least squares |29| will be used. Following the procedure of 
the bisector of the two ordinary least squares described in [29| > we nn d the best fit to be 

b = 1.7828 and A = 52.7838 , (4) 

with la uncertainties 

cr 6 = 0.0072 and er A = 0.0041 . (5) 

The best-fit calibration line Eq. with b and A in Eq. ([4]) is also plotted in the right pan el of Fig. [TJ 
From Eq. (O, it is easy to see that the calibration in this work is better than the one in [la ]. 




FIG. 1: Left panel: The Hubble diagram of 557 SNIa (black diamonds) and 50 low-redshift GRBs (red stars) 
whose distance moduli are derived by using cubic interpolation. Right panel: 50 GRBs data (red stars) in the 
log_E p ,i /(300keV) — logEiso/erg plane. The best-fit calibration line is also plotted. See the text for details. 



Next, we extend the calibrated Amati relation to high-redshift, namely z > 1.4. Since E p ^ for the 59 
GRBs at z > 1.4 have been given in [16|, [l_8|, [l9( , we can derive E- mo from the calibrated Amati relation 
Eq. @ with b and A in Eq. (jU). Then, we derive the distance moduli \i for these 59 GRBs at z > 1.4 
using Eqs. ^ and J3]) while their Sboio can be taken from [lj| [r| [l9[ . On the other hand, the propagated 
uncertainties are given by Q 

1/2 

(6) 
where 

in which <JE iso Bya is the systematic error and it accounts the extra scatter of the luminosity relation. As 
in 0, by requiring the x 2 /dof of the 50 points at z < 1.4 in the log-E Pi i /(300keV) — logi? iso /erg plane 
about the best-fit calibration line to be unity, we find that 
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cr|. = 0.1526. 

^iso.sys 



(8) 
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Note that in principle er|,. is a free parameter. However, if we permit it to vary with cosmology, 
like in e.g. [65], there might be room for the circularity problem. Even one does not care this problem, 
the constraints on cosmological models become loose, mainly due to the fact that the number of free 
parameters has been increased. On the other hand, in the cosmology-independent calibration method 
proposed in [l2j], we have not used any cosmology when we calibrate GRBs at z < 1.4, unlike in [65| 
they calibrate all GRBs without using SNIa at low-redshift. Mainly due to the nature of the calibration 
method proposed in [l2j , we have no freedom to determine a 2 E . by cosmology, and hence we must 
use the method in Q to fix it by requiring x 2 /dof = 1. We admit that this prevents us to learn the 
systematics dominating the Amati relation. Nevertheless, we plot the derived distance moduli \i with 
la uncertainties for these 59 GRBs at z > 1.4 in Fig. [21 We also present the numerical data of these 
59 GRBs in Table |TTJ For convenience, we call them "Hymnium" GRBs sample. It is worth noting that 
they are obtained in a completely cosmology-independent manner, and hence can be used to constrain 
cosmological models without the circularity problem. 




2 4 6 8 



z 

FIG. 2: The Hubble diagram of 557 Union2 SNIa (black diamonds) and 59 high-redshift Hymnium GRBs (red 
stars) whose distance moduli are derived by using the calibrated Amati relation. The dashed line indicates z — 1.4. 
See the text for details. 



III. METHODOLOGY TO CONSTRAIN COSMOLOGICAL MODELS 

Very recently, as mentioned in Sec.|H the SCP collaboration released their Union2 dataset which consists 
of 557 SNIa (20|, whereas the WMAP Collaboration also released their 7-year CMB data (WMAP7) 
in [2l|. On the other hand, as of the end of 2009, there are 109 long GRBs with measured redshift 

and spectral peak energy Ha (and [MEGS)- 

The number of usable GRBs is significantly increased. 
Motivated by these significant updates in the observations of SNIa, CMB and GRBs, it is natural to 
consider the joint constraints on cosmological models with the latest observational data. Here, we use 
the combination of 557 Union2 SNIa dataset 20], 59 calibrated Hymnium GRBs dataset (obtained in 
this work, and whose numerical data are given in Table |TTJ), the shift parameter R from the WMAP 7 
data [2l|, and the distance parameter A of the measurement of the BAO peak in the distribution of SDSS 
luminous red galaxies [22|, HH . 

The data points of the 557 Union2 SNIa compiled in [20( and the 59 Hymnium GRBs dataset in Table [TTI 
of this paper are given in terms of the distance modulus (J>obs(zi)- On the other hand, the theoretical 
distance modulus is defined as 



Hth(zi) = 51og 10 D L (zi) + i-i , 
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TABLE I: The numerical data of 50 low-redshift GRBs at z < 1.4. The first 4 columns are taken from ftfl. rt& rtgj . 
whereas the last column is derived by using cubic interpolation from the 557 Union2 SNIa. These 50 low-redshift 
GRBs can be used to calibrate the Amati relation. See the text for details. 

where /^o = 42.38 — 51og 10 h and h is the Hubble constant Hq in units of 100 km/s/Mpc, whereas 

W1 = {l+Z) [wm- (10) 

in which E = H/Hq, and p denotes the model parameters. Correspondingly, the \ 2 from the 557 Union2 
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SNIa and the 59 Hymnium GRBs is given by 



(11) 



where a is the corresponding ler error. The parameter /io is a nuisance parameter but it is independent of 
the data points. One can perform an uniform marginalization over (1q. However, there is an alternative 



way. Following [3CJ, |3l| , the minimization with respect to /io can be made by expanding the of Eq. (jTTJ) 
with respect to /io as 

where 
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Eq. (I12[) has a minimum for /io = B/C at 



xl(p) = A(p)- 
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obviously, we can instead minimize xji which is independent of fiQ. In fact, by 
this method, we have just marginalized over /io analytically, and then simplified our computing. 

There are some other observational data relevant to this work, such as the observations of CMB 
anisotropy [2l| and large-scale structure (LSS) [IH HI]- However, using the full data of CMB and LSS to 
perform a global fitting consumes a large amount of time and power. As an alternative, one can instead 
use the shift parameter R from the CMB, and the distance parameter A of the measurement of the baryon 
acoustic oscillation (BAO) peak in the distribution of SDSS luminous red galaxies. In the literature, the 
shift parameter R and the distance parameter A have been used extensively. It is argued that they are 
model-independent [32], while R and A contain the main information of the observations of CMB and 
BAO, respectively. As is well known, the shift parameter R of the CMB is defined by (321. [33l] 



R = fl 



1/2 
mO 



dz 
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(14) 



where fl m o is the present fractional density of pressureless matter; the redshift of recombination = 
1091.3 which has been updated in the WMAP7 data [21]. The shift parameter R relates the angular 
diameter distance to the last scattering surface, the comoving size of the sound horizon at and the 
angular scale of the first acoustic peak in CMB power spectrum of temperature fluctuations [H, [33|. 
The value of R has been updated to 1.725 ± 0.018 from the WMAP7 data [2l|. On the other hand, the 
distance parameter A of the measurement of the BAO peak in the distribution of SDSS luminous red 
galaxies [22| is given by 



A = tf r { 2 a E(z b ) 
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(15) 



where z b = 0.35. In |2J 
scalar spectral index n s 
the total x 2 is given by 



the value of A has been determined to be 0.469 (n s /0.98)-° 35 ± 0.017. Here the 
is taken to be 0.963, which has been updated from the WMAP7 data 21]. So, 
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(16) 



where xft is given in Eq 



M: Xcmb = (R- RobsY/'JR and Xb ao = ( A - A obe ) 2 /a A . The best-fit 
model parameters are determined by minimizing the total x 2 - As in [H, |66j], the 68% confidence level 
is determined by A% 2 = X 2 ~ Xmin — 2.3 and 3.53 for n p = 1, 2 and 3, respectively, where n p is 
the number of free model parameters. On the other hand, the 95% confidence level is determined by 



Ax 2 



X 



x 2 ■ 



< 4.0, 6.17 and 8.02 for n p = 1, 2 and 3, respectively. 
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TABLE II: The numerical data of 59 calibrated GRBs at z > 1.4. The first 4 columns are taken from [Wl. rt8l. rt§j ■ 
whereas the last column is derived by using the calibrated Amati relation. These 59 calibrated GRBs are called 
Hymnium sample, and can be used to constrain cosmological models without the circularity problem. 
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IV. OBSERVATIONAL CONSTRAINTS ON COSMOLOGICAL MODELS 



A. Parameterized models 



1. ACDM model 



As is well known, for the ACDM model, 

e{z) = ^/n m0 (i + zf + (i - n m0 ) . (17) 

It is easy to obtain the total x 2 as a function of the single model parameter fi m o for the ACDM model. 
We present the corresponding \ 2 an d likelihood C oc e~ x I 2 in Fig. [3j The best fit has Xmin = 566.173, 
whereas the best- fit parameter is fl m o — ®-2706^q'q\^1 (with la uncertainty) 1q 0263 ( w hh 2a uncertainty). 




0.15 0.2 0.25 0.3 0.35 0.4 0.15 0.2 0.25 0.3 0.35 0.4 



^m0 ^m0 



FIG. 3: The x 2 and likelihood C oc e~ x2/2 as functions of f2 m0 for the ACDM model. 



2. XCDM model 



Next, wc consider the XCDM model. It is also well known that in the spatially flat universe which 
contains pressureless matter and dark energy whose equation-of-state parameter (EoS) is a constant w x , 
the corresponding E(z) is given by 

E(z) = yfn m0 (l + z) 3 + (1 - fi m0 )(l + z) 3( - 1+w x) . (18) 

By minimizing the corresponding total x 2 in Eq. (I16|) . we find the best- fit parameters tt m o — 0.2704 and 
w x = —0.9967, while Xmin = 566.168. In Fig. |H we present the corresponding 68% and 95% confidence 
level contours in the fl m o — w x parameter space for the XCDM model. 



3. CPL model 

Here, we consider the familiar Chevallier-Polarski-Linder (CPL) model ;35], in which the EoS of dark 
energy is parameterized as 

w de = w + w a (l - a) = w Q + w a — — , (19) 

1 + 2 
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where wq and w a are constants. As is well known, the corresponding E{z) is given by [34L l36l [3 



E{z) 



n m0 (i + zf + (i - n m0 ) (i + z )3(i+»o+».) 



exp 



3w a z 
'l + z 



1/2 



(20) 



There are 3 independent parameters in this model. By minimizing the corresponding total \ 2 in Eq. (I16p , 
we find the best-fit parameters O m o = 0.2719, wo — —1.0451 and w a — 0.2635, while Xmin = 566.007. In 
Fig. [5l we present the corresponding 68% and 95% confidence level contours in the wo — w a plane for the 
CPL model. Also, the 68% and 95% confidence level contours in the f2,„o — wq plane and the £l nl Q — w a 
plane for the CPL model are shown in Fig. [51 
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FIG. 4: The 68% and 95% confidence level contours in the fi m o 
The best-fit parameters are also indicated by a black solid point. 



parameter space for the XCDM model. 



DGP model 



One of the simplest modified gravity models is the so-called Dvali-Gabadadze-Porrati (DGP) 
braneworld model [38L l39j . which entails altering the Einstein- Hilbert action by a term arising from 
large extra dimensions. For a list of references on the DGP model, see e.g. |40j,|4l| and references therein. 
As is well known, for the spatially flat DGP model (here we only consider the self-accelerating branch), 
E{z) is given by [H-l4l| 



E(z) = y/fl m0 (l + zf + fl rc + y/a r 
where ft rc is a constant. It is easy to see that E(z = 0) = 1 requires 



a 



i 



(21) 



(22) 



Therefore, the DGP model has only one independent model parameter fl rc . Notice that < f2 rc < 1/4 
is required by < fl m o < 1. It is easy to obtain the total % 2 as a function of the single model parameter 
Q rc . In Fig. [3 we plot the corresponding x 2 and likelihood C oc e~ x I 2 . The best fit has Xmin = 611.794, 
whereas the best-fit parameter is VL rc — 0.1356igQQ4g (with ler uncertainty) 1q 0099 ( w ^h 2a uncertainty). 
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FIG. 6: The same as in Fig. except for the Q m o — wo plane and the Sl m o — w a plane. 



C. New agegraphic dark energy model 

In [42|, [43[ , the so-called new agegraphic dark energy (NADE) model has been proposed recently, based 
on the Karolyhazy uncertainty relation which arises from quantum mechanics together with general 
relativity. In this model, the energy density of NADE is given by [42|, |43| 



Pq = 



3n 2 rrip 



(23) 
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FIG. 7: The y 2 and likelihood £ oc e x / 2 as functions of f2 rc for the DGP model. 



where m p is the reduced Planck mass; n is a constant of order unity; rj is the conformal time 



dt 

il = I — 

a 



da 



(24) 



in which a = (1 + z) 1 is the scale factor. Obviously, fj = 1/a, where a dot denotes the derivative with 
respect to cosmic time t. The corresponding fractional energy density of NADE reads 



o„ 



Pq 



3m 2 H 2 H 2 r) 2 



(25) 



From the Friedmann equation H 2 = (p m + p q ) / (3m 2 ) , the energy conservation equation p m +3Hp m = 0, 
and Eqs. (|2l?)) — (|2"5"|) . we find that the equation of motion for f2 g is given by [i^, l43| 



dtt q 

dz 



-0,(1-0,) 



3(l + «)- 



-VO, 



(26) 



From the energy conservation equation p q ■ 
the EoS of NADE is given by [11 EH 



3H(p q +p q ) = 0, and Eqs. ((Ml) — it is easy to find that 



3n a 



(27) 



When a — >■ oo, O g — >■ 1, thus w q — > — 1 in the late time. When a — >• 0, 9 — > 0, so 0/0 appears in w q and 
hence we cannot directly obtain w q from Eq. (|27[) . Let us consider the matter-dominated epoch, in which 



H 2 



oc a 3 . Thus, a 1 / 2 da oc dt = adr/. Then, we have rj oc a 1 / 2 . From Eq. ([23 



p 9 oc 



From 



the energy conservation equation p q + 3Hp q (l + w q ) = 0, we obtain w q = —2/3 in the matter-dominated 



Comparing w q = —2/3 with Eq. (|2T|) . 



epoch. Since p m oc a -3 and p 9 oc a" 1 , it is expected that 9 

we find that O g = n 2 a 2 /4 in the matter-dominated epoch as expected. For a <C 1, provided that n is of 
order unity, 9 <c 1 naturally follows. There are many interesting features in the NADE model and we 
refer to the original papers [42|, |43[ for more details. 

At first glance, one might consider that NADE is a two- parameter model. However, as shown in [42], 
NADE is a single-parameter model in practice, thanks to its special analytic feature 0,j = n 2 a 2 /4 = 
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n 2 (l + z)~ 2 /4 in the matter-dominated epoch, as mentioned above. If n is given, we can obtain Cl q (z) 
from Eq. (|26p with the initial condition Cl q (zi n i) — n 2 (l + Zj n j)~ /4 at any 2j n j which is deep enough 
into the matter-dominated epoch (we choose Zi n i = 2000 as in [42J), instead of £l q (z = 0) = 1 — fi TO o 
at z = 0. Then, all other physical quantities, such as fi m (z) = 1 — f2 g (;z) and w q (z) in Eq. (|27l) . can 
be obtained correspondingly. So, Q m o = £l m (z = 0), £l q o — fl q (z — 0) and w q o = w q (z = 0) are not 
independent model parameters. The only free model parameter is n. Therefore, the NADE model is a 
single- parameter model in practice. 

From the Friedmann equation H 2 = (p m + p q ) j (3m 2 ), we have 



E{z) 



fl m0 (l + zf 

i-n q {z) 



1/2 



(28) 



If the single model parameter n is given, we can obtain Vt q (z) from Eq. (|21)|) . Thus, we get f2 m o = 
1 — Q. q (z = 0). So, E(z) is at hand. Therefore, we can find the corresponding total x 2 m Eq- (HH])- In 
Fig. [51 we plot the total x 2 and likelihood C oc e~ x I 2 as functions of n. The best fit has Xmm = 594.449, 
whereas the best-fit parameter is n = 2.8865^0 0^8 (with la uncertainty) i° 1622 (with 2er uncertainty). 




FIG. 8: The X 2 and likelihood £ cx e" x2/2 as functions of n for the NADE model. 



D. Holographic dark energy model 

The so-called holographic dark energy (HDE) has been studied extensively in the literature. It is 
proposed from the holographic principle [4J, |45[ in the string theory. For a quantum gravity system, 
the local quantum field cannot contain too many degrees of freedom, otherwise the formation of black 
hole is inevitable and then the quantum field theory breaks down. In the thermodynamics of the black 
hole [4f| |47| . there is a maximum entropy in a box of size L, namely the so-called Bekenstein entropy 
bound Sbh, which scales as the area of the box ~ L 2 , rather than the volume ~ L 3 . To avoid the 
breakdown of the local quantum field theory, Cohen et al. [H| proposed a more restrictive bound, i.e., 
the energy bound. If pa is the quantum zero-point energy density caused by a short distance cut-off, 
the total energy in a box of size L cannot exceed the mass of a black hole of the same size [48| , namely 
L 3 pA < Lm 2 . The largest IR cut-off L is the one saturating the inequality. Thus, 

PA = 3c 2 m 2 p L- 2 , (29) 
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where the numerical constant 3c 2 is introduced for convenience. For the HDE model proposed in [4£ 
the cut-off L has been chosen to be the future event horizon Rh, which is given by 



f°°dt 
Rh = a J — = a 



a 



da 



(30) 



From Eqs. (j29l) . ([30]) . and the energy conservation equation p\ + 3Hpj\(l + w\) = 0, it is easy to find 
that (see e.g. (49M53j |) 



^ = -(1 + z)- l Sl A (l - SIa) ( 1 + -yfth. 
az \ c 

where £l\ is the fractional energy density of HDE. From the Friedmann equation H 2 = (p n 
we have 



E(z) 



fimo(l + Zf 



i - n A (z) 



1/2 



(31) 



Pa) I (3m 2 ) 



(32) 



There are 2 independent model parameters, namely fi m o and c. One can obtain ^Ia(z) by solving the 
differential equation (|3ip with the initial condition &a{z = 0) = 1 — O m o- Substituting Q, A (z) into 
Eq. (|32p . we can find the corresponding E(z) and then the total x 2 in Eq. (|16|l . By minimizing the total 
X 2 , we find the best-fit parameters fl m0 = 0.2764 and c = 0.7482, while xL in = 566.215. In Fig. |H we 
present the corresponding 68% and 95% confidence level contours in the £l m o ~ c parameter space for the 
HDE model. 




0.32 



FIG. 9: The 68% and 95% confidence level contours in the Q m o — c parameter space for the HDE model. The 
best-fit parameters are also indicated by a black solid point. 



E. Ricci dark energy model 



The so-called Ricci dark energy (RDE) model was proposed in [54| , which can be regarded as a variant 
of HDE mentioned above, while its corresponding cut-off L in Eq. (|2"9"| is chosen to be proportional to 
the Ricci scalar curvature radius. In [54], there is no physical motivation to this proposal for L in fact. 
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Recently, in [55[ it is found that the Jeans length Rcc which is determined by Rq C = H + 2H 2 gives the 
causal connection scale of perturbations in the flat universe. Since the Ricci scalar is also proportional 
to H + 2H 2 in the flat universe, the physical motivation for RDE has been found in [55j actually. In the 
RDE model, the corresponding p\ is given by [HJ] 



p A = 3am 2 p (H + 2H 2 



(33) 



where a is a positive constant (when L is chosen to be Rcc in Eq. (|29l) . one can see that a = c 2 in fact). 
Substituting Eq. (1331) into Friedmann equation, it is easy to find that p§i . l5~i| 



E(z) 



(l + z) 3 + 1 



(1+*) 



4,-2/ a 



-1 1/2 



(34) 



There are 2 independent model parameters, namely Q m o an d a . By minimizing the corresponding total 
X 2 in Eq. {T5J), we find the best-fit parameters il m0 = 0.3223 and a = 0.3559, while x min = 589.026. In 
Fig. I10[ we present the corresponding 68% and 95% confidence level contours in the fl m Q — a parameter 
space for the RDE model. 



0.4 



0.38 



0.36 



0.34 



0.32 




0.28 



0.3 



0.32 

^m0 



0.34 



0.36 



FIG. 10: The 68% and 95% confidence level contours in the £l m o 
best-fit parameters are also indicated by a black solid point. 



a parameter space for the RDE model. The 



COMPARISON OF MODELS 



In the previous section, we have obtained the constraints on 7 cosmological models with the latest 
observational data, namely, the combination of 557 Union2 SNIa dataset [2fJ, 59 calibrated Hymnium 
GRBs dataset (obtained in this work), the shift parameter R from the WMAP7 data [21], and the 
distance parameter A of the measurement of the BAO peak in the distribution of SDSS luminous red 
galaxies [22|, [23j. In Table UTIl we summarize the results for these 7 models. Here, we would like to briefly 
consider the comparison of these models. Naively, one can compare them with their Xmini i- e -' the best 
model has the smallest Xmin ( v i ce versa). However, as is well known, in general Xmin decreases when 
the number of free model parameters increases. So, Xmin ^ s not a g°°d criterion for model comparison. 
Instead, a conventional criterion for model comparison in the literature is Xmm/^ /' m which the degree 
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Model 


ACDM 


XCDM 


CPL 


DGP 


NADE 


HDE 


RDE 


Xmin 


566.173 


566.168 


566.007 


611.794 


594.449 


566.215 


589.026 


k 


1 


2 


3 


1 


1 


2 


2 


\min fdof 


0.918 


0.919 


0.920 


0.992 


0.963 


0.919 


0.956 


ABIC 





6.421 


12.687 


45.621 


28.276 


6.468 


29.280 


AAIC 





1.995 


3.834 


45.621 


28.276 


2.042 


24.853 


Rank 


1 


2 ~ 3 


4 


7 


5 ~ 6 


2 ~ 3 


5 ~ 6 



TABLE III: Summarizing all the 7 models considered in this work. 



of freedom dof = N — k, whereas N and k are the number of data points and the number of free model 
parameters, respectively. We present the Xmin/dof for all the 7 models in Table Mil On the other 
hand, there are other criterions for model comparison in the literature. The most sophisticated criterion 
is the Bayesian evidence (see e.g. [56j and references therein). However, the computation of Bayesian 
evidence usually consumes a large amount of time and power. As an alternative, one can consider some 
approximations of Bayesian evidence, such as the so-called Bayesian Information Criterion (BIC) and 
Akaike Information Criterion (AIC). The BIC is defined by p57j 

BIC = -2 In Umax + A; In TV , (35) 

where C max is the maximum likelihood. In the Gaussian cases, Xmin ~ ~ 21n£ rnax . So, the difference in 
BIC between two models is given by ABIC = Ax 2 rnin + AklnN. The AIC is defined by [58j 

AIC = -2ln£ max + 2k. (36) 

The difference in AIC between two models is given by AAIC = Ax 2 min + 2Ak. We refer to e.g. [5j| 
for some relevant works with BIC and AIC. In Table ITTT1 we also present the ABIC and AAIC of all the 
7 models considered in this work. Notice that ACDM has been chosen to be the fiducial model when we 
calculate ABIC and AAIC. From Table iLUl it is easy to see that the rank of models is coincident in all 
the 3 criterions (x^ n /do/, BIC and AIC). The ACDM model is still the best one, whereas DGP model 
is the worst one. The detailed rank of all the 7 models is also given in Table HTT1 

VI. CONCLUSIONS 

In the present work, by the help of the newly released Union2 compilation which consists of 557 SNIa, 
we calibrated 109 long GRBs with the well-known Amati relation, using the cosmology-independent 
calibration method proposed in [l2j]. We have obtained 59 calibrated high-redshift GRBs which can be 
used to constrain cosmological models without the circularity problem (we call them Hymnium GRBs 
sample for convenience). One can directly read off the numerical data of these 59 Hymnium GRBs from 
Table [H] of this paper. They are also available upon request to our email address. Considering their 
unique ability to extend our vision up to redshift z = 8.1, we recommend to use these 59 Hymnium 
GRBs in the relevant works. 

We considered the joint constraints on 7 cosmological models from the latest observational data, namely, 
the combination of 557 Union2 SNIa dataset [2(|, 59 calibrated Hymnium GRBs dataset (obtained in 
this work), the shift parameter R from the WMAP7 data [2l|, and the distance parameter A of the 
measurement of the BAO peak in the distribution of SDSS luminous red galaxies [22j,[23|. Note that in 
general it is better to see how the parameters are constrained both with and without GRBs. However, 
in fact the cosmological constraints on the same models without GRBs have been considered in e.g. [52j |. 
although they used slightly earlier SNIa and CMB data. So, we do not consider the constraints without 
GRBs in the present work. Comparing with the previous results in the literature (e.g. [52]), one can 
find that the constraints obtained in this work are tighter, thanks to these newly improved observational 
data. We also briefly considered the comparison of these 7 cosmological models. The ACDM model is 
still the best one from the perspective of the latest cosmological observations. 
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